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ABSTRACT 

We calculate the spectrum of photons resulting from electromagnetic cascades through 
thermal radiation, and examine the consequences of including triplet production in 
these cascades. We assume that the cascade is one-dimensional, and we find that this 
approximation is justified in the present work for thermal radiation with temperature 



less than 10 



Results are obtained for both monoenergetic and power-law primary 



spectra, and for a variety of path lengths. We find that triplet production is particularly 
important in electron-photon cascades through thermal radiation when the primary 
energy exceeds 10 5 m 2 c 4 /fcT for propagation over small path lengths. The importance 
of triplet production decreases as the path length increases, and it has no effect on 
saturated cascades. 
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^1 INTRODUCTION 

nI Electron-photon cascading through radiation fields involv- 
ONj n g photon-photon pair production (PP) and inverse Comp- 
Q ton scattering (ICS) governs the spectrum of high-energy 
O R adiation from a variety of astrophysical sources (Aharo- 
JL nian, Kirillov-Ugryumov & Vardanian 1985; Zdziarski 1986; 
^ Svensson 1987). Such cascades have already been applied 
"^in determining the spectra of 7-rays from energetic objects 
such as compact X-ray sources exhibiting very high-energy 
activity (Mastichiadis 1986; Protheroe & Stanev 1987) and 
active galactic nuclei (Mastichiadis & Protheroe 1990). They 
also occur in spatially varying anisotropic radiation fields, 
modifying the primary 7-ray spectra emitted from com- 
pact lueminous objects (Protheroe, Mastichiadis & Dermer 
1992). In addition, such cascades may enable 7-rays above 
10 15 eV to be observed from extragalactic objects at dis- 
tances far in excess of the pair production interaction length 
in the microwave background (Protheroe 1986). At even 
higher energies, cascades initiated as a result of pion pho- 
toproduction interactions of cosmic ray protons in the mi- 
crowave background may result in an observable 7-ray back- 
ground at TeV energies (Halzen et al. 1990). 
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However, at the highest energies inverse Compton scat- 
tering becomes less important than the process (7e — > 3e) 
referred to here as triplet production (TP). The importance 
of triplet production in astrophysics has been emphasized 
by Mastichiadis, Marscher & Brecher (1986), Mastichiadis 
(1991) and Dermer & Schlickeiser (1991). As was shown in 
those papers, triplet production cannot be ignored in colli- 
sions where the photon energy in the electron rest frame 
exceeds 10 2 mc 2 , and triplet energy losses dominate over 
inverse Compton energy losses for energies greater than 
10 6 mc 2 , where m is the electron rest mass. Despite this, TP 
has not yet been incorporated into full cascade calculations. 



In this paper, we consider electron-photon cascades in 
which TP acts both as a competitive energy loss mechanism 
to ICS and as a particle re-injection mechanism. In previous 
work (Mastichiadis, Szabo & Protheroe 1991) we considered 
the effects of TP on cascades through 'thin' slabs of thermal 
radiation of temperature T. We found that for ®E > 10 2 , 
where = kT and E is the primary electron energy, both 
in units of mc 2 , TP can have a significant effect on the spec- 
trum of particles emerging from the slab. Here we extend 
this work to include 'thick' slabs as well as saturated cas- 
cades, i.e. where the energies of the photons are below the 
pair production threshold on the background radiation field. 
In Section 2 we briefly describe the physical processes that 
are important in this type of cascade. The method of cal- 
culation, which employs a combination of the Monte Carlo 
technique and the matrix doubling technique of Protheroe 
& Stanev (1993), is described in Section 3. The results are 
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discussed in Section 4 and we give our conclusions in Sec- 
tion 5. 



2 THE PROCESSES 

For all three processes (a = ICS, TP, PP), the mean interac- 
tion length in an isotropic and uniform radiation field with 
temperature is given by (assuming E 3> 1) 



ra(e,0) 



<Tq(s)(l - fin) 

2 



dfide, (1) 



where n(e, 0) is the differential photon number density of 
the radiation field, cr a (s) is the cross-section for process a, s 
is the square of the total centre of momentum frame energy 
in units of m 2 c 4 , fi = cos is the cosine of the interaction 
angle, and fic is the velocity of the primary particle (J) = 1 
for PP). Here E is the energy of the primary particle and e is 
the soft photon energy, both in units of mc 2 . For ICS there 
is no threshold energy and hence eics = and /tics = 1. 
The threshold conditions for PP and TP are s > 4 and 
s > 9 respectively, and the corresponding photon thresh- 
old energies are epp = E -1 and etp = 2E _1 . Similarly, 
/tpp = 1 — 2(eE)~ 1 and /itp = 1 — 4(e_E) -1 . We compare the 
mean interaction lengths as a function of ®E for the three 
processes in Fig. 1(a). In each case, the mean interaction 
lengths have been scaled to the Thomson mean interaction 
length x T (Q) = l/[N(Q)a T ], where N(Q) = J n(e, 0)de 
is the total number density of the radiation field, and <tt 
is the Thomson cross-section. From Fig. 1(a) we can see 
that TP has a shorter mean interaction length than ICS for 
QE > 10 2 . 

The mean energy loss rates for ICS and TP are given 

by 



E 



\dt/ t a (E,Q) 



K a (E,e)Ec 
x a {E,Q) ' 



(2) 



where K a (E,Q) is the mean inelasticity for process a as a 
function of the primary electron's energy E and the temper- 
ature of the radiation field 0. We have calculated K a (E,Q) 
for both ICS and TP using a Monte Carlo calculation, which 
will be described in Section 3.1. In the case of TP where 
there are three electrons (throughout this paper 'electrons' 
is taken to mean both electrons and positrons) in the final 
state, we have taken the total energy of the two lowest en- 
ergy electrons to be the energy lost. The resulting energy 
loss time-scales, t a (E, 0), are shown in Fig. 1(b) in units of 
the time-scale for Thomson scattering. For comparison, we 
have also added the interaction rate for PP to Fig. 1(b). We 
can see that, while TP interactions occur more frequently 
than ICS interactions for QE > 10 2 , TP energy losses only 
dominate over ICS energy losses for QE > 10 6 . This is be- 
cause the mean inelasticity for TP is relatively small (Mas- 
tichiadis 1991) and falls off as E~° 65 for QE > 10 2 , whereas 
for ICS the mean inelasticity is approximately constant over 
the same energy range. 



3 METHOD OF CALCULATION 

The application of Monte Carlo techniques to cascades dom- 
inated by the physical processes described above presents 
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Figure 1. (a) The logarithm of the mean interaction length 
[scaled by a:x(©)] f° r ICS (solid curve), PP (dotted curve) and 
TP (dashed curve), as a function of log(G£/). (b) The reciprocal 
energy loss time-scale for ICS and TP [scaled by c/x'y^Q)], as a 
function of ®E. Also shown is the interaction rate for PP. Curves 
are defined as in part (a). 



some difficulties, which we will try to address in the follow- 
ing sections. In Section 2 we saw that, while TP occurs more 
frequently than ICS for QE > 10 2 , relatively small amounts 
of energy are lost per interaction. Hence high-energy elec- 
trons will undergo many triplet interactions, in which very 
little energy is lost, until they undergo a catastrophic ICS. 
If one follows all the produced pairs, the computing time 
required is increased considerably. 

The approach we use here is based on the matrix mul- 
tiplication method described by Protheroe (1986) and sub- 
sequently developed by Protheroe & Stanev (1993). We use 
a Monte Carlo program to calculate the yields of secondary 
particles due to interactions with the thermal radiation. The 
yields are then used to build up transfer matrices which de- 
scribe the spectra of particles produced after propagating 
through a slab of radiation with thickness Ax. Manipulation 
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of the transfer matrices as described below then enables one 
to calculate the spectra of particles resulting from propaga- 
tion through slabs of arbitrarily large thickness. 

3.1 Particle yields 

We define the quantity Y^ 13 as the probability of producing 
a particle of type fi =e (electron) or 7 (photon) in the bin 
centred on energy Ej, when a primary particle with energy 
in the bin centred on E t undergoes an interaction of type a 
with the radiation field. To calculate Y" 13 we use a Monte 
Carlo simulation. For ICS and PP we have used the com- 
puter code described by Protheroe (1986, 1990), updated 
to model interactions with a thermal photon distribution 
of arbitrary temperature. For TP, we have used the code 
described by Mastichiadis (1991) to model the interaction, 
and we have sampled the initial conditions of the interaction 
as described by Mastichiadis et al. (1991). While calculat- 
ing the yields for ICS and TP, we also calculated the mean 
inelasticity for ICS and TP used in Section 2 by summing 
over the fractional energy lost by the primary electron dur- 
ing each interaction. 

3.2 Transfer matrix calculation 

The transfer matrix, T^" (Arc), gives the number of secon- 
daries of type v =e or 7, with energy in the bin centred on 
Ej, emerging from a slab of thermal photons of thickness Ax 
when a primary of type fi =e or 7 and energy in the bin cen- 
tred on Ei is injected into the slab. To calculate the transfer 
matrices we have used a modification of the semi-analytical 
technique described by Protheroe & Stanev (1993). If Ax 
is much shorter than the shortest interaction length in the 
cascade, i.e. Ax <C it, then 
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where l/x e (Ei) = l/xics(E t ) + 1/xtp(-E 8 ) gives the total 
mean interaction length for e7 interactions, and xics(Ei), 
xpp(_E 8 ) and xtp(-E 8 ) are given by equation (1). 

3.3 Matrix doubling method 

The method used to calculate the spectrum of particles 
emerging from a slab of arbitrary thickness is described by 
Protheroe & Stanev (1993), and is summarized below. Once 
the transfer matrices have been calculated for a slab of thick- 
ness Ax, the transfer matrix for a slab of thickness 2 Ax is 
simply given by applying the transfer matrices twice, i.e. 
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k = j 



T7=(Ax)T e J(Ax) 



+ T77(Ax)T77(Ax)]. (10) 

The new matrices may then be used to calculate the transfer 
matrices for a slab of thickness 4 Ax, and so on. A slab of 
thickness 2" Ax only requires the application of this 'matrix 
doubling' n times. The spectrum of electrons and photons 
that escape from the slab is then given by 

CO 

F !(*) = ^[n e (x)^(0)+T7;(x)f7(0)], (11) 

i = 3 

CO 

where ff(0) and F?(0) are the input electron and photon 
spectra (number of electrons or photons per unit energy in 
the bin centred on Ei), and x = 2" Ax. In this way, cascades 
through thick slabs can be modelled quickly and efficiently. 

The matrix method and matrix doubling technique have 
been used for many years in radiative transfer problems (van 
de Hulst 1963; Hovenier 1971), and in such problems scatter- 
ing of the radiation means that 'reflection' in each layer Ax 
must also be properly taken into account. Reflection is ne- 
glected in the present work because the angular deflection of 
the particles produced in the three interactions considered is 
negligible, and even the cumulative deflection after repeated 
interactions in the whole cascade is usually small, as we will 
show below. 

3.4 Angular deflection of cascade particles 

In this section, we justify neglecting the angular deflection 
of particles in interactions, and show that it is unnecessary 
to include reflection terms in the matrix calculations. We do 
not give here a rigorous treatment of the angular distribu- 
tion of particles emerging from an interaction, as this is quite 
complicated. Instead, we make order of magnitude estimates 
of the maximum possible deflections. In practice, the deflec- 
tions will usually be less than we estimate below. Generally, 
we will proceed along the following lines. We estimate the 
scattering angle per interaction, Aip, and then estimate the 
typical number of interactions required to 'cool' the particle, 
N . An estimate of the cumulative deflection after cascading 
is then ip ~ N~ Aip. Note that the number of scatterings, N, 
will in general be much smaller than the slab thickness over 
which we are calculating the cascade divided by the mean 
interaction length, except for very small slab thicknesses for 
which deflection will in any case be unimportant. 
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We start with a consideration of electrons with energies 
below the TP threshold. These electrons will cool exclusively 
by ICS, and the interactions will be in the Thomson regime. 
For scattering of a soft photon of energy e by a relativis- 
tic electron of energy E, the typical energy in the electron 
rest frame before scattering is e ~ Ee. Since scattering is in 
the Thomson regime, this also equals the typical energy after 
scattering. The scattering is isotropic (within a factor of 2) in 
the electron rest frame, and so the transverse component of 
the photon's momentum after scattering is of order e . This 
will also be the order of magnitude of the transverse com- 
ponent of the electron's momentum after scattering, in both 
the electron rest frame and the laboratory (LAB) frame. 
Since the electron's energy loss per interaction is small when 
in the Thomson regime, the angular deflection per interac- 
tion is of order Aip ~ e' / E ~ e. In the present case of 
blackbody radiation of temperature 0, we have e ~ 2.70, 
and hence Aip ~ 2.70. The typical energy loss of the elec- 
tron per interaction is E 2 e. The number of interactions to 
'cool' the electron is then typically N ~ E / E 2 e ~ 1/2. 70£. 
Thus the cumulative scattering angle is 

ip ~ N^Aip ~ (13) 
V&E 

As expected, this increases as the energy decreases. At the 
lowest energy used in the present work, for which ®E = 
10 -4 , the cumulative scattering angle will be small (less than 
10°) provided that < 10 -3 (i.e. for temperatures less than 
6 x 10 6 K). In reality, the limit will be somewhat higher 
because our 'typical' deflections are in fact overestimates. 

For PP well above threshold (i.e. s ^> 4), the deflection 
angles are quite small as we will show below. By integrating 
over the differential cross-section for PP (Jauch & Rohrlich 
1976), we find that the median production angles in the 
centre of momentum (CM) frame of the two electrons are 
Aib* ~ 2/s* and (ir — 2/s*) with respect to the initial di- 
rection of the energetic photon in the CM frame. The angle 
between the energetic photon's direction in the CM frame 
and its direction in the LAB frame is of order y/s/E and so 
can be neglected in comparison with the median production 
angles. The transverse component of a produced electron's 
momentum is p± ~ (y/s/2) sin Aip* ~ s". The longitudi- 
nal component of momentum in the LAB frame is obtained 
by a Lorentz transformation using 7cm — E/y/s. We obtain 
jo 1 1 ~ E(l — 1/y/s) (corresponding to the high-energy particle 

with Aib* ~ 2/s*) and py ~ E / y/s (corresponding to the 
low-energy particle). We see that usually the high-energy 
particle carries off most of the energy. The LAB frame pro- 
duction angles are then Aip ~ s< / E (high-energy particle) 
and Aib ~ s* / E (low-energy particle). 

For ICS in the extreme Klein-Nishina regime, by inte- 
grating over the differential cross-section (Jauch & Rohrlich 
1976) we find the median angle in the electron rest frame be- 
tween the photon's direction after scattering and its initial 
direction to be Aip' ~ y/2/( eE )i . The angle between the 
photon's direction in the electron rest frame before scatter- 
ing and the boost direction is ~ ir — 1/E, and this may be ne- 
glected in comparison with the median scattering angle. The 
median LAB-frame angle between the scattered photon's di- 
rection and the initial direction of the electron is found by 
Lorentz transformation to be Atb ~ y/2(eE)i/E. With the 



help of further Lorentz transformations one obtains the me- 
dian LAB-frame energies after scattering of the electron, 
Ex ~ E/VeE, and of the photon, e t ~ E(l - 1/y/lE), 
and the median LAB-frame electron scattering angle Aib ~ 
V2(eE)i/E. Since, for ICS, s = 1 + 2Ee(l - fi), we find 
that the median LAB-frame angles and energies of the pho- 
ton and electron after scattering are within a factor of 2 of 
those of the high-energy and low-energy electrons in pair 
production interactions well above threshold. 

For ICS in the extreme Klein-Nishina regime, it is found 
that the high-energy particle taking most of the energy is al- 
most always the scattered photon, and so in the cascade well 
above the PP threshold the high-energy particle alternates 
between an electron and a photon. The high-energy particle 
has negligible deflection in comparison with the low-energy 
particle. In the cascade, particles will migrate from high to 
low-energy by a variety of routes - in some collisions a par- 
ticle will come out as the low energy particle, while in other 
it will come out as the high-energy particle. First we con- 
sider the cooling of the 'high-energy particle', i.e. we follow 
only the highest energy particle coming out of each interac- 
tion. The energy lost per interaction is AE ~ E/y/s, giv- 
ing N ~ y/s collisions to cool the particle, with a resulting 
cumulative deflection of ib ~ y/s/E. If instead we follow 
the low-energy particle produced, as it is effectively cooled 
in one interaction we have ib ~ s*/E. This latter deflec- 
tion imposes the stronger constraint and we see that, if the 
small-angle criterion is satisfied at some energy, it will auto- 
matically be satisfied at higher energies. We will postpone a 
discussion of the implications of this result until we discuss 
PP near threshold. 

In the case of PP (and also TP) the largest deflection 
occurs close to threshold when the particles are produced 
almost isotropically in the CM frame. We shall consider 
the 'worst case scenario' of the produced particles sharing 
equally the energy available, and being produced in the CM 
frame with their velocity vectors transverse to the primary 
particle's direction. In this case, the energy of an electron 
in the CM frame is y/s/n, where n is 2 (3 for TP). Ne- 
glecting the electron mass, which will again cause overesti- 
mation of scattering angles, the transverse momentum (in 
any frame) is also y/s/n. The longitudinal component of the 
LAB-frame momentum is ~/cMy/s/n, giving a deflection an- 
gle of Aib ~ y/s/E. Interactions such as these occurring near 
threshold will result in particles produced below the thresh- 
old for further cascading-type interactions (PP, or ICS in 
the Klein-Nishina regime), and in any case the mean in- 
teraction lengths near threshold are extremely large and 
so such interactions are rare. Hence N ~ 1 and we obtain 
ib ~ y/s/E. At threshold y/s = 2, and therefore for small- 
angle deflections (ib < 10°) we require E > 12. For PP, how- 
ever, y/s = y/2Ee(l — fi) ~ 3y/@E with head-on collisions 
(p ~ —1) more likely at threshold, and so we obtain the con- 
straint < 0.03 on the temperature of the radiation in or- 
der that angular deflections be small. With this temperature 
constraint, the small-angle criterion is clearly also satisfied 
in the high-energy part of the cascade, as discussed above 
both for the case of the deflection of the 'high-energy par- 
ticle' and for individual scatterings when a particle emerges 
as the low-energy particle. 

Finally, we consider TP well above threshold. An up- 
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per limit to the deflection per interaction is given by Aip < 
AE/E, where AE is the average energy lost per interaction. 
The number of interactions required to coo l an electron is 
typically N < E/AE. Hence ij> < ^/AE/E = ^/xtpAtpc. 
Taking the interaction length and energy loss time-scale for 
TP from Fig. f, we see that, at energies where the TP in- 
teractions occur more frequently than IC (@E > fO 2 ), ip is 
small. 

We conclude that angular deflection is definitely unim- 
portant in the present work provided that the temperature 
of the radiation is < I0 -3 , and is probably unimportant 
to somewhat higher temperatures. As the temperature ap- 
proaches ~ I, however, pair production by photon-photon 
collisions of photons of the radiation field, as well as inter- 
actions in any associated plasma, becomes important, and 
our results will not be valid. The neglect of reflection terms 
due to scattering in our matrix calculation is justified when 
< I0 -3 . However, for higher temperatures (but still with 
<C 1 ), our results presented below in terms of the spectra 
of particles or radiation emerging 'from a slab of thickness 
x' should instead be interpreted to mean spectra of parti- 
cles or radiation 'at a time x/c after injection of the primary 
particle into the radiation field'. 



4 RESULTS AND DISCUSSION 

We have calculated the photon spectrum emerging from 
slabs of various thicknesses for both monoenergetic pri- 
maries and a power-law primary spectrum. The power- 
law spectrum has the form Fo(QE) = A(QE)~ T , where 
-Emin < E < -Emax and A is fixed by normalizing the spec- 
trum to contain one particle. Results are given below for 
both electron and photon primaries. 



4.1 Monoenergetic injection 

The spectra of electrons and photons that escape from a rel- 
atively thin slab of thickness 2 5 xt(0) are shown in Fig. 2. In 
(a) primary electrons of energy E = IO 5 /0 are injected into 
the slab, and we can see that the inclusion of TP causes 
a significant difference in both the shape and the normal- 
ization of the resultant spectra. Both the photon and the 
electron spectra are enhanced by as much as four orders 
of magnitude at energies below E ~ IO 3 /0. Similarly in 
Fig. 2(b) spectra are shown for primary electrons of energy 
E = fO 7 /©- H ere the effect of including TP is even more 
prominent, with the spectra being enhanced by as much as 
eight orders of magnitude at low energies. 

In Fig. 3 we investigate the effect of increasing the 
slab thickness on the escaping photon spectrum. We show 
in Fig. 3(a) the photon spectra obtained for electron pri- 
maries with an energy of E = IO 5 /0 after cascading through 
slabs of thickness 2 1o xt(0) (thinnest curves), 2 15 xt(0) and 
2 2 % T (0) (thickest curves). In each case the results with TP 
included (solid curves) are compared with those where TP is 
not included (dotted curves). From Fig. 3(a) we can see that, 
while the slab is 'thick' (2 10 xt and 2 15 xt for this primary 
energy), the inclusion of TP can cause the normalization 
of the photon spectrum below the PP threshold to be en- 
hanced by up to several orders of magnitude. However, once 
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Figure 2. The spectrum of particles and radiation emerging 
from a slab of thickness 2 5 X'y(&) f° r monoenergetic primaries. 
Note that, because of the way in which E and G are defined, 
(@E) 2 F(@E) is dimensionless. Results are shown for the case 
where TP is included (solid curves) and the case where it is not 
(dotted curves). The thick (thin) curves represent the spectra of 
photons (electrons and positrons) that escape from the the slab, 
(a) Primary electron of energy E = 10 5 /B. (b) Primary elect ron 
of energy E = 10 7 /B. 



the cascade approaches saturation (e.g. 2 20 xt for this pri- 
mary energy) the results for the two cases become identical 
below the PP threshold. Similar behaviour can be seen in 
Fig. 3(b) where we show the photon spectra emerging from 
slabs of thickness 2 15 xt(0) (thinnest curves), 2 2 % T (0) and 
2 25 x T (0) (thickest curves) for primary electrons of energy 
E = IO r /0. Here the inclusion of TP may enhance the pho- 
ton spectrum at low energies by as much as four orders of 
magnitude. In each case, the cascade is seen to saturate for 
smaller slab thicknesses when TP is included, and there is 
also a significant difference in the shape of the spectrum be- 
tween the PP threshold and the primary electron's energy. 

We have also calculated the photon spectrum produced 
when photon primaries instead of electrons are injected into 
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Figure 3. The spectrum of radiation emerging from slabs of var- 
ious thicknesses for monoenergetic primaries. Results are shown 
for the case where TP is included (solid curves) and the case 
where it is not (dotted curves), (a) Primary electron of energy 
E = 10 5 /G for slabs of thickness 2 10 x^^Q) (thinnest curves), 
2 15 xt(@) and 2 20 x^(Q) (thickest curves), (b) Primary electron 
of energy E = 10 7 /G for slabs of thickness 2 15 x T (@) (thinnest 
curves ),2 20 x T (&) and 2 25 :E T (e) (thickest curves), (c) As for part 
(b), but with photon primaries. 



the slab. Results for slabs of thickness 2 15 xt(0) (thinnest 
curves), 2 2O xt(0) and 2 25 xt(0) (thickest curves) and E = 
lO r /0 are shown in Fig. 3(c). Again TP is seen to have 
an important effect for cascades through 'thick' slabs. The 
magnitude of the enhancement below the PP threshold is 
less than for the equivalent case using electron primaries 
because the 7-rays pair-produce, and hence the mean energy 
of the first electrons and positrons produced is less than 
E = lO 7 /0- 

The above results can be understood qualitatively as 
follows. Any electron in the system with E® ^> 10 2 will 
interact with the ambient photons predominantly by TP in- 
teractions. As was emphasized by Mastichiadis (1991), in 
the absence of ordinary PP these interactions will tend to 
steepen the resulting photon spectrum and enhance the nor- 
malization over the 'pure' ICS case. The inclusion of PP 
should not change the above picture as long as the slab is 
not too thick. For this type of slab, 7-rays have a finite prob- 
ability of escaping. However, those 7-rays that interact pro- 
duce an electron-positron pair which will produce a number 
of secondary pairs through TP interactions. The secondaries 
are produced with rather low energies, and as such they will 
cool by ICS, producing additional photons below the PP 
threshold. As a result, the final photon spectrum will have 
a normalization that will be higher than that expected in 
the usual ICS-PP cascades. Furthermore, one would expect 
the ratio of the two normalizations to be given by the pair 
multiplicity for TP (the number of secondaries produced per 
primary electron) as defined by Mastichiadis (1991). There- 
fore, by increasing the primary energy, one would expect 
the difference from the usual cascades to be more signifi- 
cant. The Monte Carlo results given above seem to support 
this notion. 

The situation changes when we consider saturated cas- 
cades. In this case, the inclusion of TP should not make any 
difference as the particular characteristics of the cascade are 
totally washed away. Both pair-producing processes, despite 
their different properties, will effectively degrade all high- 
energy particles and photons to energies below threshold. 
Therefore the end result should be the same irrespective of 
the way in which the particles lose their energy inside the 
slab. 

It is evident that in going from moderately thin slabs, 
where the inclusion of TP can produce a dramatic differ- 
ence of several orders of magnitude over the usual cascades, 
to very thick slabs, where the two cases are practically iden- 
tical, one has a range of slab thicknesses where the inclusion 
of TP is of decreasing importance. 



4.2 Power-law injection spectrum 

The results for a power-law primary electron spectrum 
with T = 1.5, £ min = lO _4 /0 and £ max = 10 7 /© are 
shown in Fig. 4(a). Escaping photon spectra are shown for 
slabs of thickness 2 15 xt(0) (thinnest curves), 2 2 % T (0) and 
2 25 x T (0) (thickest curves). Again an enhancement in the 
low-energy photon spectrum is seen, although in this case 
the inclusion of TP only increases the normalization by a 
factor of ~ 3 for 2 15 xt and 2 20 xt- As for the case of mo- 
noenergetic injection, once the cascade becomes saturated 
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Figure 4. The spectra of radiation emerging from slabs of var- 
ious thicknesses for a power-law primary spectrum of the form 
F (@E) = A(@E)- r , £ min < E < S max . Results are shown 
for the case where TP is included (solid curves) and the case 
where it is not (dotted curves). In each case the input spectrum 
is indicated by the dashed curve, (a) Primary electron spectrum 
with T = 1.5, £ min = 10~ 4 /e and £ max = 10 7 /© for slabs of 
thickness 2 15 ^ T (B) (thinnest curves), 2 20 x T (&) and 2 25 ^ T (B) 
(thickest curves), (b) As for part (a), except with T = 2.0. (c) As 
for part (a), but with photon primaries. 



(~ 2 25 xt in this case) the low-energy photon spectra be- 
come identical. We have also calculated results for a steeper 
primary electron spectrum (T = 2.0), which are shown in 
Fig. 4(b). In this case, the inclusion of TP has little effect 
on the results, as most of the low-energy photon spectrum 
is dominated by cascades induced by primary electrons with 
energy below E ~ 10 6 /© where ICS dominates the energy 
losses. 

The escaping photon spectra that result when a power- 
law photon spectrum (T = 1.5, -E m in and _E m ax as before) 
is injected into the slab are shown in Fig. 4(c). While an 
enhancement below the PP threshold is observed, its mag- 
nitude is again less than for the equivalent case with electron 
primaries. 



5 CONCLUSIONS 

In the present paper we have examined the effect of triplet 
production (TP) on electron-photon cascades through a 
thermal radiation field. These cascades have until now been 
examined by assuming that the only relevant processes are 
inverse Compton scattering (ICS) and photon-photon pair 
production (PP). By applying a Monte Carlo code and sub- 
sequently the matrix doubling technique, we have shown 
that the inclusion of TP can have significant effects on the 
emerging photon spectra. The importance of TP becomes 
more pronounced for high primary energies and for relatively 
thin slabs. The differences in the spectra occur both below 
and above the PP threshold. As the slab thickness increases 
the differences below the PP threshold diminish, and for the 
limiting case of an infinite slab (i.e. a saturated cascade) the 
inclusion of TP will have no practical consequences. How- 
ever, we find that the inclusion of TP causes saturation to 
occur for smaller slab thicknesses, and so even for ultra-thick 
slabs, where the spectra are almost identical below the PP 
threshold, there will still be significant differences above the 
PP threshold. The significance of TP in astrophysical envi- 
ronments will thus depend solely on the individual source 
parameters, which in turn determine whether a saturated 
cascade occurs. 
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